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Abstract. We analyze the colour magnitude diagram of the field population in the direction of the bulge cluster 
NGC 6553 to discuss the stellar field component projected towards the cluster. In particular we investigate the 
nature of a peculiar feature, called "red flare" which appears as a faint turnoff, located 2 magnitudes about below 
the turnoff of the cluster, with a sparse trace of a sub-giant branch just above. A detailed modeling shows that 
the red flare is consistent with the expected location of the turnoff and sub-giant branch of the bulge population, 
once that the appropriate metallicity range and spatial distribution of the bulge stars are adopted. From the data, 
hints can be derived about the spatial shape and the basic stellar parameters of the bulge. 
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1. Introduction 

The structure of the Galactic Bulge and the nature of its 
stellar population have been controversial problems in the 
last decade. Most of the current knowledge is based on 
models derived from low resolution, photometric infrared 
observations. The data interpretation requires some hy- 
potheses on the distribution of the dust and some assump- 
tion on the disk/bulge separation. Efforts have been as 
well devoted to the interpretation of field color-magnitude 
diagrams (CMDs) projected in the direction of the inner 
bulge along the minor and the major bulge axis. These dia- 
grams appear quite complex, with a nu mber of unexpected 
featur es. Re cently Feltzing & Gilmore (2000) and Bcaulicu 
et al. (^00 iD published HST CMDs of the globular cluster 
NGC 6553, projected at galactic coordinates I — 5 and 
b = —3. In their diagrams they detect a peculiar feature 
called "red flare" . This feature appears as a faint turnoff, 
located 2 magnitudes about below the turnoff of the clus- 
ter, with a sparse trace of a sub-giant branch just above. 
The two papers leave the problem open, suggesting dif- 
ferent interpretations. Feltzing & Gilmore (2000) suggest 
a possible contamination by the background Sagittarius 
galaxy. Beaulieu et al. (2001) seem more inclined to a 
bulge background turnoff. They fail to reproduce that 
turnoff with an isochrone of the same [Fe/H] content than 
the cluster, namely [Fe/H]= —0.4 concluding than the 
mean bulge population must be more metal poor. Both the 
interpretations, however, have to face some difficulties: the 



Sagittarius galaxy hypothesis implies a wide contamina- 
tion, not present in the other fields, while the "more metal 
poor" bulge field one is based on data where cluster and 
field stars are mixed together, preventing any firm con- 
clusion. In addition Beaulieu et al. (2001) reproduce the 
bulge with a single metallicity stellar population located at 
8 Kpc distance from the Sun, while recent papers based 
on high resolution spectroscopy suggest a wide range of 
[Fe/H], from —1.6 to 0.55 with a mean value of —0.14 
(Rich & Mc William 2000). As far as the cluster is con- 
cerned, the most recent literature determinations, based 
on high resolution spectroscopy, give [Fe/H]= —0.55 ±0.2 

Fe/H]= -0.16 ± 0.1 
revised upward to 



or Z 
or Z 
Fc/H ] 



0.01 (Barbuy |200C| ), and 
0.02 (Cohen et al 
-0.06 ± 0.15 or Z 



1999) 
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0.027 by Carretta et al. 
( 200l| ) with evidence of an exces s of a — process elements. 
From these results Cohen et al. ( 1999| ) and Carretta et al. 
( ^00 1| ) conclude that the value of [Fe/H] for NGC 6553 
is comparable to the mean of the Galactic Bulge. This 
implies that there is no strong evidence in the recent lit- 
erature supporting an higher metallicity of the cluster in 
comparison to the field. For this reason we decide to pro- 
ceed to a deeper investigation aimed to detect the true na- 
ture of this feature using the galactic models developed in 
Padova (Vallenari et al. 200C) and the most rece nt CMDs 
of NGC 6553 published by Zoccali et al. ( |200l| ) , where 
relative proper motions are provided. This allow us to dis- 
card cluster objects from field stars. The goal is to discuss 
whether an appropriate admixture of bulge and disk stars 
can be responsible of this red flare. 
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In Section 2 the Galaxy model is presented, in Section 
3 the observed and simulated CMDs are shown and the 
resulting shapes of the bulge are discussed. Finally the 
conclusions are drawn in Section 4. 

2. Modeling the Galaxy 

The simulation of the Gala xy is d one with the code already 
descri bed by Bertelli et al. ( 1995 ) and revised by Vallenari 
et al. ( 20001 ). The generation of the synthetic population 
makes use of the set of stellar tracks by Salasnich et al. 
(2000) for isochrones with a enhanced abundances. We 
adopt for the old disk population an age of 10 Gyr while 
the youngest population is reaching an age of 2 Gyr. 

Evidence have been advanced in literature that the 
bulge might not have spherical symmetry. Several mod- 
els have been proposed, going from spherical to triaxial 
symmetry. In our galaxy model the following choices are 
included: 

1) C OBE-DIRBE oblate spheroidal model (Dwek et 
al. |1995| , their model GO). It is a Gaussian type function, 
with axial ratios 1:1:0.56, similar to the values derived on 



the basis of the near infrared maps by Kent et al. (1991). 

2) COB E-DIRBE constrained triax ial m odel by Kent 
et al. (1991) as modified by Dwek et al. (1995) (their model 
G2). It is a "boxy" Gaussian models, with axial ratios 
1:0.22:0.16. 



3) the triaxial model introduced by Binney et al. ( 1997 ) 
to fit COBE-DIRB E dat a with the best fitting parameters 
by Bissantz et al. ( 1997 ) (thereafter BO). The axial ratios 
are 1:0.6:0.4. 



4) an exponential triaxial function (Dwek et al. 1995 
model E2) based on the distribution of the IRAS Mira 
variable. The axial ratios are 1:0.18:0.39. This model is 



suggested by Dwek et al. ( 1995 ) as giving one of the most 
rehable fit to the COBE-DIRBE observations. 

The angle between the sun-centre line and the major 
axis of the bar is ranging from 25° (Bissantz et al. 1997) 
to 20° (Dwek et al. |1995| G2 model) while it becomes 36° 
in the case of the E2 model. 

A stellar population having ages from 12 to 10 Gyr and 
metal content ranging from Z = 0.008 to 0.05 is adopted 
as suggested by the most recent papers based on high 



resolution spectroscopy (Rich & McWilliam 2000, Barbuy 
et al. p^99|). 



3. The red flare population in the field of NGC 
6553 

The first step towards the interpretation of the data is 
the separation of the field from the cluster stars. Zoccali 



et al. (2001) compare images of NGC 6553 taken at dif- 
ferent epochs in order to derive star proper motions. As 
expected, these authors find that the cluster and the 
field populations follow two different distributions of rel- 
ative proper motions and can consequently be disentan- 
gled. Cluster objects are distributed by construction in a 
Gaussian centered on relative proper motion in I and h 



(dl, dh) = (0, 0) while field population is clustered around 
(dZ, dh) = (-0.2 45, -0 .017) pixels. Following the analysis 
by Zoccah et al. ( ^OOl h, in their CMD of the region around 
NGC 6553 (see Fig.|l) we select stars having dl < -0.2 
(corresponding to a proper motion of fii < 5.89 mas yr~^ 
of the cluster with respect to the field objects). Those 
stars are believed to sample the foreground and back- 
ground population. Consequently, the field population is 
disentangled from the cluster objects. The corresponding 
CMD is presented in Fig.|^. About 1200 stars are selected 
in the field. In particular the red flare is brought into ev- 
idence. Since NGC 6553 is located in the direction of the 
Galactic bulge, the most obvious hypothesis is that this 
population is simply due to a superposition along the line 
of sight of mainly bulge stars with a residual disk contam- 
ination. The stars redder than {V — I) ^ 2.15 would then 
be sub-giant and clump objects of disk and bulge. In or- 
der to reproduce the red flare under this hypothesis, it is 
necessary to remember that a single isochrone at a fixed 
distance of 8 Kpc is not an appropriate representation of 
the bulge population. In fact the bulge is an admixture of 
objects of various metallicity, and has a spatial distribu- 
tion as discussed in previous Sections. Fig. || presents the 
predicted distributions of the bulge stars with the distance 
from the Sun in the line of sight of NGC 6553 at varying 
bulge models. The vast majority of the models presented 
in literature expect bulge stars to be located at a distance 
from 5 to 11 Kpc from the Sun. Fig. ^ where isochrones 
are superposed on red flare data shows how populations 
of various metallicity placed at different distance from the 
Sun can reproduce the location and the colour of the sub- 
giant branch and clump. 

A Galactic disk component is needed to reproduce the 
main sequence region bluer than (V — I) 1.8 and roughly 
brighter than V = 20.5 mag, as far as the appropriate red- 
dening is adopted (see the following Section) . The evolved 
stars of this disk population are located in the same CMD 
region of the sub-giant branch and clump stars of the 
bulge. 



3.1. Extinction determination 



It has been proved (Vallenari et al. 2000 ) that the blue 
edge of the main sequence in the CMD of the disk and 
bulge population can be interpreted in terms of extinction 
along the line of sight. Comparing the observed main se- 
quence edge location with theoretical simulations at vary- 
ing extinction, we derive the values of Ay along the line of 
sight as a function of the heliocentric distance (see Fig.^). 
Using this method, we find that at the distance of 5 Kpc 
from the Sun the absorption is Ay = 2.7 mag. This value 
is larger than Ay — 2.08 obtained by Zoccali et al. (2001) 
from the fit of the cluster NGC 6553, located at 5.3 Kpc 
from us. However, these authors claim they analyze a low 
reddening region inside the cluster. In fact, differential red- 
dening is present in the field. This is evident from Fig.^ 
where cluster data of the whole field obtained by Zoccali 
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Fig. 1. The observational CM d iagram of the cluster NGC 
6553 from Zoccali et al. (2001). Superposed on it a grid 
of isochrones at varying reddening from Ay — 2.0 to 2.8 
in steps of 0.1. The age is 12 Gyr, the metal content Z — 
0.019, with enhanced abun dance s of a elements in respect 
to Solar (Salasnich et al. 2000 ). The assumed distance 
modu lus is (to — M)o = 13.6 as derived by Zoccali et al. 
2001. 



et al. (2001) are plotted together with isochrones at vary- 
ing extinction. If the width of the cluster sub-giant branch 
is mainly due to differential reddening, the data are con- 
sistent with Av ranging from 2.1 to 2.7. This result is in 
substantial agreement with the analysis made by Sagar et 
al. (|199|), Guarnieri et al. (pJ97|), Beaulieu et al. (pOOl|) 



who find a differential reddening of about 0.5 mag inside 
the field. As far as the foreground and background pop- 
ulation of Fig.|^ is concerned, the main contribution to 
the absorption is due to the disk interstellar gas and dust 
at distances closer than 5-6 Kpc. In fact, to reproduce the 
colour of the main sequence field population at higher dis- 
tances, an additional extinction of 0.1 is required leading 
to a final value of Ay = 2.8. 

3.2. The results 

Using the model described in Section ^ including thin 
disk and bulge we reproduce the CMD of Fig.|^. The re- 
sult of the simulation is shown in Fig.^. It is quite obvi- 
ous looking at Fig.^ that different spatial distributions of 
the stars in the bulge can predict different turnoff magni- 
tudes. The termination point magnitude changes at max- 
imum of 0.3-0.4 mag because of this effect. A detailed 
discussion can be found in Vallenari et al. (1999). Here 
we recall that "rounder" distributions (model GO or BO) 
result in brighter magnitudes (~ 20.4), while flat distri- 



butions having large position angle (model E2) produce 
fainter turnoff (^ 20.7). The luminosity function of the 
disk population reaches a maximum at V 20.5 mag cor- 
responding to a distance from us of about 7 kpc. Because 
of the combined effect of disk and bulge population, to 
compare observational and simulated CMDs is then not 
trivial. We make use of a test (see Vallenari et al. 2000 



for more detail on the adopted procedure), is defined 
as = ^/i^J2ii^i ~ /o-i where and bi are the ob- 
served and the expected number of stars per magnitude 
bin, and v is the number of degrees of freedom (which is 
equal to the number of bins minus 1, since the in the sim- 
ulations we impose that the total number of model stars 
is equal to the observed total number of stars). In con- 
sequence of this definition, high probability models have 
~ 1. Due to the poor statistics, all the models, but the 
model E2 are in good agreement with the data, resulting 
in a ~ 1-5 — 2. Although it is difficult to reach a firm 
conclusion, models BO and GO seem to result in a more 
convincing agreement with the data, model G2 is not in 
disagreement, while the model E2 is definitely less prob- 
able (x^ '^3). As a final comment on the adopted bulge 
population, we point out that while the metallicity and 
distance spread are fundamental to reproduce the obser- 
vational width of the sub-giant branch, the age range is 
poorly constrained by the data, due to the poor statis- 
tics. This analysis strongly points in favor of the fact that 
the red fiare is simply due to the bulge and disk popu- 
lations. From our simulations, we estimate that the disk 
contributes at maximum for the 30% to the red sub-giant 
branch and clump population. 

4. Conclusions 

Using simulations derived from models of the Galaxy and 
the recent HST CMDs of NGC 6553 and its projected field, 
we are able to reproduce the bulge and disk field contri- 
butions and the so-called "red flare" discussed in Feltzing 
& Gilmorc ( |200C| ) and Beauheu et al. ( |200l| ). This feature 
turns out to be a natural consequence of bulge and disk 
field contamination. The spread in distance and metallic- 
ity is fundamental to reproduce the observational CMD 
of the field. The bulge is modeled using an age from 10 
to 12 Gyr and a metallicity ranging between Z — 0.008 
and Z = 0.05. The cluster is fitted with an isochrone of 
Z — 0.019 which is inside the metallicity range of the 
bulge field. We recall however that while constraints set 
on metallicity are stringent, the age range is poorly deter- 
mined, due to the limited statistics. 
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Fig. 2. The observational CM diagram of the field popu- 
lation towards NGC 6553 together with three isochrones 
of Z = 0.019, age 12 Gyr, and distance of 5, 8,10 Kpc 
(solid lines), and an isochrone of Z = 0.008, age 12 Gyr 
and distance of 5 Kpc (dotted line) 
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Fig. 3. The expected distribution of bulge stars with the 
distance from the Sun at varying spatial distributions in 
the line of sight of NGC 6 553. O pen squares indicate the 
model by Bissantz et al. ( 1997 ), circles show the model 
E2, crosses refer to the model G2 and finally filled squares 
represent the model GO (see text for detail). 



Fig. 4. The simulated CM diagram of the background and 
foreground population towards NGC 6553. This simula- 
tion includes bulge and disk stars. 
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